Spectral properties of spherically confined dusty plasma crystals 
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A combined theoretical and experimental analysis of the normal modes of three-dimensional sphe- 
rially confined Yukawa clusters is presented. Particular attention is paid to the breathing mode and 
the existence of multiple monopole oscillations in Yukawa systems. Finally, the influence of dissipa- 
tion on the mode spectrum is investigated. 

PACS numbers: 52.27.Lw, 52.27.Gr 



I. INTRODUCTION 



Spherically confined dust crystals have recently attracted large attention, see e.g. [T] and |2] for an overview. After 
the exploration of the ground state structure [31 |U |5J [B] and of metastable states and their probability 13 |H1 H] a 
question of particular interest is the dynamical excitation spectrum. The normal modes are a key property describing 
the response of a finite system to external excitation and also the melting behavior, e.g. [ini E]. Particularly 
interesting is the situation in dusty plasmas where the normal mode spectrum of two-dimensional crystals could be 
analyzed experimentally, e.g. [T^lIIS]- Here, we extend this analysis to three-dimensional spherical Yukwa crystals 
(Yukawa balls) and present experimental results. From the theory side, we recall the determination of the normal 
modes of confined systems. We concentrate on one of the key modes - the breathing mode (BM) describing uniform 
expansion and contraction of the whole system. Recently it was shown [14j that a "true" (i.e. uniform) BM does 
not exist in Yukawa systems. Instead, there may exist several similar modes which, however, all deviate either from 
uniform or radial motion of the particles. The brief derivation of Ref. |14j is reproduced here in a more detailed way. 
One of the important consequences - the possible existence of multiple monopole modes - is demonstrated. Finally, 
we take into account the effect of dissipation on the normal mode spectrum which is crucial for dusty plasmas. 



II. NORMAL MODES OF FINITE CLUSTERS 



The crystals are characterized by a 3-dimensional, classical system of N identical particles harmonically confined 
by the potential 4>{r) — miJ^r^ (2 and interacting with a Yukawa-type pair interaction v(r) — exp(— Kr)/r. The 
hamiltonian is then given by (^' indicates no summation over equal indices) 



_m 

i—1 i—1 i,j — l 



t/(rGR3«) 

To investigate the spectral properties of dusty plasma crystals at low temperatures we consider small excitations from 
a ground state or metastable state r* = (r^, r2 . . . r^) e M.^^ , which is a minimum of the potential energy and thus 
fulfills the equations 

N 

= V.t/(r)|_^ = + ^<^V. < N. (2) 

The small excitation is a time-dependent function r{t) in the configuration space with |r(t) — r*| <C 1. Due to the 
small amplitude of the oscillation the potential can be approximated harmonically by 

U{r)KU{r*) + ]^{r-r*fH'''{r-r*) (3) 
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where we used ^ and the definition of the Hessian matrix — W"^C/(r)|r=r-*- Because is a real symmetric 
3iV X SA'^-matrix and positive semidefinite as well, its eigenvalue problem |14j 

Xmr = H'''r (4) 

defines 3 A'' eigenvalues Xj > and 3 A'' linearly independent eigenvectors rj, which form a basis in the configuration 
space and may be conveniently chosen orthonormal. Within this basis the excitation can be expanded 

3N 

r{t)^r* +J2cj{t)r,, (5) 

so that the time-dependence is fully determined by the coefficients Cj(t) - the normal coordinates. These normal 
coordinates obey equations of motion, which follow from the hamiltonian ([ij 

0^mr + VC/(r) i mr + Ti''* (r - r*) 



3N 3N 3N (^Q^ 

^mY,Cj (t) r, + (i)^'* ^1 ^ E [^J- (^) + ^J^J (*)] ■ 



Because the eigenvectors are linearly independent each normal coordinate Cj(t) fulfills — Cj(t) + XjCj(t) with the 
solution 



Cj (t) = Aj cos ( ^/Xj t + Bj )yj < 3N, (7) 

in which the constants Aj and Bj have to be determined from the initial conditions r(0),f(0) of the excitation. 
Thus in general, the excitation relative to the ground state, r{t) — r*, is a superposition of oscillating motions 
Aj cos{y^t + Bj)rj, called normal modes. Each such normal mode describes a collective motion of all particles 

oscillating with the same frequency ujj — ^/Xj - the normal frequency. The respective displacements of the particles 
are given by the eigenvector r, . Consequently, in order to investigate the spectral properties of the dusty plasma 
crystals we have to investigate the normal modes, and hence have to consider the eigenvalue problem Q of the 
Hessian H.^ . 

For harmonically confined Coulomb systems (k — 0) detailed theoretical studies have been performed for d— 1, 2, 3 
dimensions, see [101 E] and references therein. It was shown that there exist three (partially degenerate) normal 
modes, which are independent of the particle number N: 

1. The d{d — l)/2 rotational modes with A = 0, which correspond to a rotation of the whole system and reflect the 
axial symmetries of the confinement. 

2. The d center of mass oscillation modes with X — uj^ expressing that the center of mass motion is independent 
of the interparticle forces. 

3. The breathing mode (BM) with A = 3cjg, which describes a uniform radial expansion and contraction of all 
particles. 

The existence of these three modes is illustrated for the two-dimensional system with = 3, 4, 5 particles in figure [T| 
where all modes of these systems corresponding to the ground state configuration are shown. 



III. NORMAL MODES IN YUKAWA SYSTEMS 



For harmonically confined Yukawa systems the interaction potential and thus the Hessian TV depend on a screening 
parameter n. Hence, in general, the normal modes will depend on this parameter as well. Within figure |2]this situation 
is shown. There, the frequencies of all normal modes are shown for the ground state of harmonically confined Yukawa 
systems with iV = 16 particles and a screening parameter ranging from K, — QioK — 20djr^ [21^. Moreover, the figure 
shows, that there are two mode frequencies independent from screening - corresponding to the rotational modes and 
the center of mass oscillation modes. These modes exist also in a Yukawa system (and have K-independent frequencies) 
due to the symmetries of the systems [15]. However, in general this is not the case for the breathing mode, which 
has no K-independent frequency. Actually, the general existence of this mode is a special property of Coulomb-like 
interactions not featured in Yukawa systems [T^ . 
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FIG. 1: (Color online) All normal modes of two-dimensional harmonically confined Coulomb systems with A'' = 3, 4, 5 particles. 
The dots picture the particles within a ground state configuration, and the arrows show the direction of the oscillatory motion. 
The Af-independent modes, i.e., the rotational modes, the modes of center of mass oscillation (c.o.m.), and the breathing modes 
are highlighted. 




FIG. 2: (Color online) K-dependence of the normal modes of a harmonically confined Yukawa system with A = 16 particles vs. 
screening parameter. The rotational modes, the center of mass oscillations and the breathing mode are highlighted dash-dotted 
(red), dashed (green), and (blue) dot, respectively. 



The existence condition of the BM becomes clear by analyzing the eigenvalue problem (|4| . The standard procedure 
to evaluate this equation is to find the roots of its characteristic polynomial yielding all eigenvalues and -vectors and to 
examine whether there exists an eigenvector t^m corresponding to a BM, i.e., describing a uniform and radial motion 
r^M oc r* . However, the degree of this equation is d-N which prohibits an analytical calculation of the eigenvalues 
and -vectors and, consequently, prohibits general statements about the existence of the BM. Therefore, we apply a 
different approach ~ the direct mode approach ~ which focuses on just one mode: within the eigenvalue equation we 
directly use the eigenvector of the breathing mode, which thus leads to the existence condition of this mode. 

Additionally, we evaluate Q within its particle components i < N, because the potentials and v are given for 
these, and obtain 



N N 

1=1 1=1 



(8) 
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Using the isotropy of 4> and the distance dependence of v one obtains for each component of the Hessian 

ViVjU{r)\r=r' - 

>'(|rrl). 



-X3 + 



„*|3 



(|rr|0"(|rr|)-<^'(|rr|)) 
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(9) 



k=l 



wherein X3 denotes the three-dimensional identity matrix. Thus ^ reduces to 

n|r*|)„ 



Amr, = 
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i\r*u\) 



(10) 



Now we use the aforementioned direct mode approach, i.e., we set the eigenvector r equal to the unit eigenvector of 
the BM tbm = \r*\~'^r* , which consequently leads to the existence conditions of this mode 



A 



BM^r* = (t>"{\r*\)r* 



N 

1=1 



4l)4Vz<iv. 



(11) 



These conditions show |14j . that a breathing mode only exists in harmonically confined systems with interaction 
potentials of the form (7 G M^o) or log(r), or in some highly symmetric systems. 

The breathing mode and its existence is closely related to the monopole oscillation (MO) [TT], which in strongly 
correlated systems is associated with the oscillation of the mean square radius R^{t) = X]t=i ^iW^ — '''{^Y 
pUl [T^ [T5] . This MO is particulary important since it can be easily excited selectively by variation of the trap 
frequency wq, e.g., by a rapid change of the frequency wo of the equilibrated system corresponding to the excitation 
given by r(0) = (1 + 0^* and f (0) = 0. 

In the case that a system has a BM, the MO describes the same motion as the BM and thus has the same 
frequency. However, in the opposite case, as for most Yukawa systems, there are several monopole oscillations with 
different frequencies. This property is directly apparent in the spectrum <I>^2(w) of the mean square radius motion 



3N 



2N- 



3N 



Cj{t)Cjit), 



(12) 



which may be derived by the square of its Fourier transform <i>^2(^) = l-^/j'2(w)p. For the mentioned excitation the 



normal coordinates are given by Cj(t > 0) 
(^ ^ 1) the spectrum is given by 



^r* ■ r, cos(wjt) and Cj(t < 0) = ^r* • r,-. Thus, for small amplitudes 



2e_ ^ {r,-r*fu:',{r,-r*f^l 

j,k=l V 



(13) 



If the BM exists, one eigenvector rj is given by tbmi which is parallel to r*, whereas all other eigenvectors are 
orthogonal to r* . Hence the spectrum contains only a single peak at the breathing frequency. In contrast, if no BM 
exists, generally many eigenvectors Vj will have a non- vanishing projection on r* [cf. figurejsjb)], so that the spectrum 
contains several MO with different frequencies wmo = and amplitudes. 

Within figure [s] a) such a spectrum is shown for a Yukawa ball with = 16 particles and a screening strength 
of K = 4.42 (ic~^. One clearly sees two peaks, thus this system contains two different monopole oscillations. The 
corresponding frequencies are = 2.34905cl;o and ~ 2.35954ti;o- 



IV. EXPERIMENTAL RESULTS 



A very intriguing experimental realization of harmonically confined (screened) Coulomb systems are the so-called 
Yukawa balls in dusty plasmas. There, a finite number of plastic microspheres are dropped into the plasma of a 
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FIG. 3: a) Part of the spectrum {^) containing two monopole modes, obtained from numerical NMA for a Yukawa ball 
with A'' = 16 particles and a screening strength of k = 4.42 dc~^. The dashed lines indicate the position of the poles of equation 
(13 1, b) In case of a vanishing BM no eigenvector is parallel to r* , but generally many eigenvectors, e.g., rj,rt,ri, have a 
non-vanishing projection on it giving rise to a (differnt) monopole oscillation. 



parallel plate rf discharge. In the plasma, the particles arrange in concentric spherical shells (see e.g. [71[H1[T2]) which 
is expected to be the ground state configuration of harmonically trapped 3D charged-particle systems [3j IH [6] . 

The harmonic confinement in the plasma is realized by a combination of various plasma forces on the microspheres. 
In the vertical direction, the action of gravity is compensated by an upward force due to the electric field in the space- 
charge sheath above the lower electrode and an upward thermophoretic force due to heating of the lower electrode. 
Horizontally, the particles are confined by a glass cube placed onto the electrode. This results in a 3D harmonic 
confinement of the charged microspheres (see e.g. for details). The charge on the typically used microspheres 
of 3.5 lira diameter is of the order of 10^ elementary charges, the interparticle distance is of about 500 /im, and the 
shielding length in the plasma is of the same order as the interparticle distance. Thus, screening strengths of k ~ 1 
are typical in these systems |S][71|S] justifying the term "Yukawa balls" for our systems. 

In the experiment, the normal modes of Yukawa balls are derived from the thermal Brownian motion of the particles 
around their equilibrium positions. The thermal particle motion is recorded with a stereoscopic video camera setup 
and the full 3D dynamics is reconstructed [H I18j . Experimentally, the normal mode spectra are obtained as the 
spectral power density Si{lu) from the Fourier transform for the mode number £ = I . . . 3N by |13j 



T 



dt 



(14) 



Here, Ci{t) = r{t) ■ ri are the projections of the Brownian particle fluctuations with velocity r{t) onto the eigenvector 
of mode £, i.e., the normal mode oscillation pattern of this mode (see e.g. Ref. ^]). 

The experimental mode spectrum is shown below for a Yukawa ball with = 10 particles. In equilibrium, the 
particles arrange on a single spherical shell as shown in Fig. [4] As an example, also the breathing mode obtained 
from the experimental particle positions is presented. The 3D coherent radial particle motion is clearly seen (In our 
Yukawa balls, the deviation from a pure Coulomb interaction is not that strong, thus the breathing mode is not much 
distorted here). 

The experimental normal mode spectrum Si{ljj) is shown in Fig. [5j It is seen that the mode frequencies are 
concentrated in a narrow band at very low frequencies of about 1-2 Hz. The dynamics of Yukawa balls occurs, 
due to the high mass of the microspheres, on very long time scales enabling the observation with video cameras. 
Also, the mode spectrum closely follows the expected mode frequencies from the mode analysis. The theoretical 
mode frequencies effectively depend only on the strength of the confinement in Eq. ([ij. By fitting the expected 
eigenfrequencies to the observed spectral power density a nice agreement is obtained and the confinement frequency 
is obtained. From the physical dimensions of the Yukawa ball and luq the dust charge is determined as Z w 900 for 
our 3.47 /im particles [19]. 
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FIG. 4: a) Bond structure model of a 3D dust cluster (Yukawa ball) with A*' = 10 particles. Note the clear spherical structure, 
b) Breathing mode for this Yukawa ball derived from the experiment. 
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FIG. 5: Experimental normal mode spectrum for a Yukawa ball with A*' — 10. Here, the spectral power density is shown in 
gray scale (dark colors correspond to high power densities). The white dots indicate the best fit of the theoretical normal mode 
frequencies to the spectral power density. 



V. NORMAL MODES IN THE PRESENCE OF DISSIPATION 



The above presented theoretical investigations of normal modes neglect several effects present in the experiment - 
most of all finite-temperature and dissipation. While inclusion of the former exceeds plain mechanics and requires 
statistics the latter can be taken into account by extending the equations of motion (|6| by a friction force —mvir with 
J' > 0, and we obtain 

ZN 

= mr + mvir + VU{r) = m'^[cj{t) + ucj{t) + XjCj{t)] Vj. (15) 
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FIG. 6: (color online) a) Illustration of the damped dynamics of a normal mode Cj{t) for three values of v with the condition 
Cj{t < 0) = 0. b) The corresponding spectra The vertical dashed line indicates the peak position in the undamped 

case, c) Damping dependence of all eigenfrequencies (i^max denotes the peak position of the spectrum) of a Yukawa Ball with 
N = 5 and a = ldc~^ ■ Darker lines highlight degenerate normal modes. 



Hence, each normal coordinate fulfills = 'cj{t) + vcjit) + XjCj{t) analogous to the single damped harmonic oscillator, 
and now depends on the damping strength i^. The solution reads {ujj = ^/Xj) 



cj{t) = A, 




for V < 2ujj 
for J/ = 2ujj 
for V > 2aj,- 



(16) 



Again, the constants Aj and Bj have to be determined from the initial conditions r(0), f(0) of the excitation. The 
three cases correspond to normal modes, which are underdamped, critically damped, and overdamped, respectively. 
Within figure [6] a) the dynamics of one mode j are illustrated for these three cases for the condition rj{0) = 0. 

The modifications due to damping can nicely be seen within the spectrum of one normal mode, which is given in 
case of the aforementioned excitation by 



27r(z.2^2 + (^2_^2)2) 



(17) 



In the undamped case {v = 0) this spectrum just contains a single distinct peak centered at the normal mode frequency. 
In contrast, dissipation causes a red-shift as well as a broadening of this peak. These effects are shown in figure [6] b). 
The frequency shift of the peak at Wmax in dependence on damping is displayed in figure |6] c) in more detail for the 
example of a Yukawa Ball with = 5 particles. 

Finally, in Fig.[7]we present the complete power spectrum for different damping parameters, for the example N — 10. 
Again, the red shift and peak broadening are evident. At sufficiently large damping more and more normal modes 
vanish, and an increasing fraction of the power density is concentrated around lo = 0. These theoretical results show a 
very similar behavior as the experimental data of Fig. |5] and allow for a sensitive diagnostics of relevant experimental 
parameters, including the particle charge, cf. Sec. |IV[ A detailed theory-experiment comparison will be presented in 
a forthcoming paper. 



VI. SUMMARY AND OUTLOOK 



Within this work we presented a theoretical and experimental analysis of the normal modes of three-dimensional 
spherially confined Yukawa clusters. Based on the predictions of Ref. |14j . we analyzed the breathing mode in more 
detail, showed the close relation to the existence of multiple monopole oscillations in Yukawa systems and presented 
new numerical results. Finally, we investigated the influence of dissipation on the mode spectrum. 
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FIG. 7: (color online) Theoretical normal mode spectrum \n\ for the cluster TV = 10 and ndc = 1 for four values of the 
damping u/ujo indicated in the figure. As initial conditions we chose Cj{t = 0) = Co for all modes. For comparison, the 
undamped spectum is shown in all figures by black dots. 
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